Primary y-ray spectra in 44 Ti of astrophysical interest 
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Primary y-ray spectra for a wide excitation-energy range have been extracted for 44 Ti from particle-y co- 
incidence data of the 46 Ti(p,/y) 44 Ti reaction. These spectra reveal information on the y-decay pattern of the 
nucleus, and may be used to extract the level density and radiative strength function applying the Oslo method. 

Models of the level density and radiative strength function are used as input for cross-section calculations 
of the 40 Ca(a, y) 44 Ti reaction. Acceptable models should reproduce data on the 40 Ca(a, y^Ti reaction cross 
section as well as the measured primary y-ray spectra. This is only achieved when a coherent normalization 
of the slope of the level density and radiative strength function is performed. Thus, the overall shape of the 
experimental primary y-ray spectra puts a constraint on the input models for the rate calculations. 

PACS numbers: 21.10.Ma, 25.20.Lj, 27.40,+z, 25.40.Hs 



I. INTRODUCTION 

The titanium isotope 44 Ti is of great astrophysical interest, 
since it is believed to be produced in the inner regions of core- 
collapse supemovae and in the normal freeze-out of Si burning 
layers of thermonuclear supernovae \ l\ , with a large variation 
of yields depending on their type [2 1. The determination of the 
44 Ti yield might reveal information on the complex explosion 
conditions. The production yield of 44 Ti directly determines 
the abundance of the stable 44 Ca, and also influences the 48 Ti 
abundance through the feeding of 48 Cr on the a chain. The 
theoretical prediction of the 44 Ti production in core-collapse 
supemovae is sensitive to the chosen reaction network and the 
adopted nuclear reaction rates. 

Cassiopeia A (Cas A) is the youngest known Galactic su- 
pernova remnant and is, at present, the only one from which 
y-rays from the 44 Ti decay chain ( 44 Ti — ► 44 Sc — > 44 Ca) have 
been unambiguously detected [2 1. The discovery of the 1 157- 
keV y-ray line from 44 Ca was reported by Iyudin et al. 13), 
while measurements of the 67.9-keV and 78.4-keV lines from 
44 Sc were presented by Renaud et al. (2). From the combined 
y-ray flux, the half-life of 44 Ti, and the distance and age of the 
remnant, an initial synthesized 44 Ti mass of i.6±gj x lfr 4 
Mq was deduced. This is thought to be unusually large, a fac- 
tor of 2 — 10 more than what is typically obtained by current 
models (see, e.g., fl4jE|). 

The main production reaction of 44 Ti is the 40 Ca(a, y) 44 Ti 
reaction channel with a g-value of 5.127 MeV; the cross sec- 
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tion for this reaction is very important to estimate the Ti 
yield. Recent cross-section measurements on this reaction |6| 
have lead to an increase of the associated astrophysical rate, 
giving a factor of ~ 2 more in the predicted yield of 44 Ti. 
Thus, the theoretical models become more compatible with 
the Cas A data; however, this would make the problem of 
"young, missing, and hidden" Galactic supernova remnants 
even more serious: supernovae that should have occured after 
Cas A are still not detected by means of y-ray emission from 
the 44 Ti decay chain. It is still an open question whether the 
Cas A is a peculiar case (asymmetric and/or a relatively more 
energetic explosion), or that the Cas A yield is in fact "nor- 
mal", since it is in better agreement with the solar 44 Ca/ 56 Fe 
ratio (2) . 

There are many uncertainties connected to the yield es- 
timate, the most severe ones being due to astrophysical is- 
sues. However, there is also significant uncertainties related 
to the nuclear physics input, in particular the nuclear level 
density (NLD), the radiative strength function (RSF) and the 
a-particle optical model potential. All these quantities are en- 
tering the calculation of the astrophysical reaction rates. The 
NLD is defined as the number of nuclear energy levels per en- 
ergy unit at a specific excitation energy, while the RSF gives 
a measure of the average (reduced) transition probability for 
a given y-ray energy. Both quantities are related to the aver- 
age decay probability for an ensemble of levels, and are indis- 
pensable in a variety of applications (reaction rate calculations 
relevant for astrophysics or transmutation of nuclear waste) as 
well as for studying various nuclear properties in the quasi- 
continuum region. 

In this work, we have extracted primary y-ray spectra from 
the decay cascades in 44 Ti populated via the two-neutron pick- 
up reaction 46 Ti(p,fy) 44 Ti. These spectra are measured for 
a wide range of initial excitation energies below the neutron 
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threshold, and they put a constraint on the functional form of 
the NLD and RSF. 

In Sec. [Il] we will describe the experimental details, give 
an overview of the analysis, and present the obtained primary 
y-ray spectra. Different normalizations of the NLD and RSF 
are discussed in Sec. [TIT] and the applied models are tested 
against the primary y-ray spectra in Sec. IV In Sec. [V] var- 
ious input models are used for estimating the 40 Ca(a, y) 44 Ti 
cross-section and reaction rates, and the calculated results are 
compared to existing data from Nassar et al. [6| and Vocken- 
huber et al. (T). Finally, a summary and concluding remarks 



are given in Sec. VI 



II. EXPERIMENT AND ANALYSIS 

The experiment was conducted at the Oslo Cyclotron Lab- 
oratory (OCL), where the Scanditronix cyclotron delivered a 
32-MeV proton beam bombarding a self-supporting target of 
46 Ti with mass thickness 3.0 mg/cm 2 . The beam current was 
« 0.5 nA and the experiment was run for ten days. Unfor- 
tunately, the target was enriched only to 86.0% in 46 Ti. The 
main impurities were 48 Ti (10.6%), 47 Ti (1.6%), 50 Ti (1.0%), 
and 49 Ti (0.8%). The reaction of interest, 46 TiQ?,fy) 44 Ti, has 
a Q-value of —14.236 MeV [7|. Particle- y coincidences from 
this reaction were measured with eight collimated Si AE — E 
particle detectors and the CACTUS multidetector system lf8l . 

The Si detectors were placed in a circle in forward direc- 
tion, 45° relative to the beam axis. The front (AE) and end 
(E) detectors had a thickness of w 140 jim and 1500 jim, 
respectively. The CACTUS array consists of 28 collimated 
5" x 5" Nal(Tl) crystals for detecting y-rays. The total effi- 
ciency of CACTUS is 15.2(1)% at E y = 1332.5 keV. Also a 
Ge detector was placed in the CACTUS frame to monitor the 
experiment. The charged ejectiles and the y-rays were mea- 
sured in coincidence event-by-event. 

To identify the charged ejectiles of the reactions, the well- 
known AE — E technique is used. In Fig. [T] the energy de- 
posited in the AE detector versus the energy deposited in the 
E detector is shown for one particle telescope. Each "banana" 
in the figure corresponds to a specific particle species as indi- 
cated in the figure. 

By gating on the triton "banana", the events with the re- 
action 46 TiQ?,fy) 44 Ti were isolated. The 46 Ti(/?,t/y) 45 Ti data 
have been published previously in Ref. [9|. The singles and 
coincidence triton spectra are shown in Fig. [2] We note that 
the excited 1.904-MeV 0+ state (E, ss 16 MeV) is rather 
weakly populated. This is expected on the basis that I — 
states are in general much weaker populated in the (p,t) reac- 
tion at 45° compared to, e.g., 1 — 2 and 4, see Ref. ifTol . The 
populated spin range is estimated to J w — 6/z based on the 
observed triton peaks and their coincident y-decay cascades, 
in accordance with the findings in Ref. iflOl . 

Two-neutron pick-up on the 48 50 Ti impurities in the target 
give rise to peaks at higher triton energies than the ground 
state of 44 Ti due to their lower reaction thresholds; the Q- 
values are -12.025 and -10.560 MeV for the 44 Ti(p, t y) 46 Ti 
and 50 Ti(p,fy) 48 Ti reactions, respectively. This means that 
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FIG. 1: (Color online). Identification of particle species using the 
AE — E technique. 
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FIG. 2: (Color online). Singles (blue) and coincidence (red) triton 
spectra. The impurities of 48 50 Ti are easily seen. The resolution is 
£3 270 - 330 keV. 



there is a background from such events in the spectra that can- 
not be removed (an w 14% effect). 

Using reaction kinematics and the known (2-value for the 
reaction, the measured triton energy was transformed into ex- 
citation energy of the residual nucleus. Thus, the y-ray spectra 
are tagged with a specific initial excitation energy in 44 Ti. Fur- 
ther, the y-ray spectra were corrected for the known response 
functions of the CACTUS array following the procedure de- 
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FIG. 3: (Color online). The first-generation spectra for each 
excitation-energy bin in 44 Ti. The dashed lines are limits set for the 
further analysis (see text). 



scribed in Ref. [ll j. The main advantage of this correction 
method is that the experimental statistical uncertainties are 
preserved, without introducing any new, artificial fluctuations. 

Information on NLD and RSF can be extracted from the 
distribution of primary y-rays, that is, the y-rays that are 
emitted first in each decay cascade. To separate these first- 
generation y transitions from the second and higher-order gen- 
erations, an iterative subtraction technique is used lfl2l . The 
y-ray spectra /, for excitation-energy bin i obviously contain 
all generations of y-rays from all possible cascades decaying 
from the excited states of this bin. The subtraction technique 
is based on the assumption that the spectra /)<, for all the en- 
ergy bins Ej < Ej contain the same y-transitions as /, except 
the first y-rays emitted, since they will bring the nucleus from 
the states in bin i to underlying states in the energy bins j. This 
is true if the main assumption of this technique holds: that the 
decay routes are the same whether they were initiated directly 
by the nuclear reaction or by y decay from higher-lying states. 

The obtained first-generation matrix P(E,Ey) of 44 Ti is 
shown in Fig. [3] The diagonal where E = Ey is clearly seen. 
Here, two peaks are particularly pronounced: one from the 
decay of the first excited state at E = Ey = 1083 keV, and the 
other from the decay of the second 2 + state at 2887 keV to the 
ground state. There are also more diagonals visible, e.g., for 
Ey = E— 1083 keV where the decay goes directly to the first 
excited 2 + state, and so on. 

Before extracting the NLD and RSF from the P matrix, 
we have set a lower limit for the y-ray energies (£™ n ), and 
a lower and upper limit for the excitation energy (£mi n , £ maK ). 
The limits on the excitation-energy side are put to ensure that 
the spectra are dominated by decay from compound states 
(£mi n ), and that the statistics is not too low {E mwi ). On the 
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FIG. 4: Original (a), unfolded (b) and first-generation (c) y-ray spec- 
trum of 44 Ti for excitation energy £ = 6.1 MeV. 



y-ray energy side the limit is set to exclude possible left-overs 
of higher-generation decay, that might not be correctly sub- 
tracted in the first-generation method (see Refs. lfl"3l[T4l and 
references therein for more details). In Fig. |4] the original, 
unfolded, and first-generation y-ray spectrum of 44 Ti for exci- 
tation energy E = 6.1 MeV are displayed. 

With the first-generation matrix properly prepared, an iter- 
ative procedure is applied to extract NLD and RSF from the 
primary y-ray spectra. The method is based on the assumption 
that the reaction leaves the product nucleus in a compound 
state, which then subsequently decays in a manner that is in- 
dependent on the way it was formed, i.e. a statistical decay 
process lfl31 . This is a reasonable assumption for states in the 
quasicontinuum, where the typical lifetime of the states is of 
the order of 10~ 15 s, whereas the reaction time is 10~ 18 s. 
In addition, the configuration mixing of the levels is expected 
to be significant if the level spacing is comparable to the resid- 
ual interaction 1131 . This is normally fulfilled in the region of 
high level density. 

If compound states are indeed populated, the y decay from 
these states should be independent of the reaction used to 
reach them. In previous works (e.g. Ref. [16|) it has been 
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demonstrated that the direct reactions ( 3 He, 3 He') and ( 3 He,oc) 
into the same final nucleus do produce very similar decay 
cascades. Also, the extracted NLD and RSF was found to 
be equal within the expected fluctuations 1 This indicates that 
even though a direct reaction such as (p,t) is used, compound 
states are likely to be populated at sufficiently high excitation 
energy (as mentioned already for the E m i n limit set in the first- 
generation matrix P(E,Ey)). 

The ansatz for the iterative method is fl3l : 



P(E,E r )~p(Ef)^(E r ), 



(1) 



meaning that the first-generation matrix P[E,Ey) is assumed 
to be separable into two vectors that give directly the func- 
tional form of the level density at the final excitation energy 

E f = E-E 

given Ey. This is done by minimizing 



y, and the y-ray transmission coefficient & for a 



1 



gmax 

I 



E 

£ 



P th (E,Ey)-P{E,Ey) 
AP(E,Ey) 



(2) 



where A^ ree is the number of degrees of freedom, AP(E,Ey) 
is the uncertainty in the experimental first-generation y-ray 
matrix P(E,Ey), and the theoretical first-generation matrix is 
given by 



P t h(E,E r ) 



P{E-Ey)J r {Ey) 
=£ min P(E — Ey) S? (Ey) 



(3) 



Every point of the p and & functions is assumed to be an 
independent variable, so that the reduced % 2 of Eq. |2} is min- 
imized for every argument E — Ey and Ey. Note that is 
independent of excitation energy according to the Brink hy- 
pothesis ifTTl . If we had an excitation-energy dependent y-ray 
transmission coefficient, 2F = 3?(E,Ey), it would in principle 
be impossible to disentangle the level density and the y-ray 
transmission coefficient. 

It is well known that the Brink hypothesis is violated when 
high temperatures and/or spins are involved in the nuclear re- 
actions, as shown for giant dipole resonance (GDR) excita- 
tions in Ref. [18] and references therein. However, since both 
the temperature reached (T « \/Ef) and the spins populated 
(/ ~ — 6fi) are rather low for the experiment in this work, 
these dependencies are assumed to be of minor importance in 
the excitation-energy region of interest here. 

To inspect how well the iterative procedure works, we have 
compared the experimental first-generation spectra for several 
excitation energies with the ones obtained by multiplying the 
extracted p and ST functions. The result for a selection of 
primary y-ray spectra is shown in Fig. [5] As can be seen, the 
agreement between the calculated and the experimental first- 
generation spectra are in general quite good, although there 
are local variations where the calculated spectra are not within 
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the error bars of the experimental ones. These variations could 
well be due to large Porter- Thomas fluctuations [19], as there 
are relatively few levels in this nucleus. 

The iterative procedure to obtain the level density and the 
y-ray transmission coefficient uniquely determines the func- 
tional form of p and £T\ however, identical fits to the experi- 
mental data is achieved with the transformations ITJl 



p{E-Ey) = Aexp[a(£ -E r )]p(E -E 
#(E r ) = Bexp{aEy)£?(E r ). 



(4) 
(5) 



Thus, to obtain the absolute normalization of the level density 
and y-ray transmission coefficient, the transformation param- 
eters A, a, and B must be determined independently. 



III. NORMALIZATIONS AND MODELS FOR 
LEVEL DENSITY AND RADIATIVE STRENGTH 
FUNCTION 

A. Level density 

Usually, the level density p(Ef) is normalized to known, 
discrete levels at low excitation energy, and to the total level 
density at the neutron separation energy p(S n ), which is de- 
duced from neutron resonance spacings Do and/or D\ (see, 
e.g., Ref. 11201 ). However, for the 44 Ti case, no neutron (or 
proton) resonance data are known since the target nuclei for 
the neutron and proton capture reactions are unstable. We are 
therefore left with only two types of experimental constraints, 
namely the known levels at low excitation energy and the ob- 
served first-generation spectra (which also depend on the y- 
ray strength). 

We have chosen two different approaches to normalize our 
data, and, out of those, to estimate the sensitivity on our re- 
sults to this normalization procedure: (f) apply theoretical 
level densities based on microscopic calculations, and (if) use 
a standard closed-form formula with global parameters. For 
the first approach, we have used recent calculations of Goriely, 
Hilaire and Koning ED (hereafter labeled GHK). For the sec- 
ond approach, we have applied the constant-temperature (CT) 
formula 112211: 



Pct(£) = j, exp I - 



(6) 



where the nuclear temperature T = 1 .50 MeV and the energy 
shift Eq = —0.08 MeV are taken from the global parameteri- 
zation of Refs. [23, 24]. The level-density data normalized to 
these two approaches are shown in Fig. [6] 

We observe that our data follow closely the known, discrete 
levels [25] at low excitation energy, and especially in the re- 
gion 2.4 < E < 5.2 MeV. This is gratifying, as it implies that 
the Oslo method does indeed give reasonable results for the 
level density. We also see that there is a decrease in the level 
density for 4.0 < E < 5.2 MeV and an abrupt increase for 
5.8 < E < 6.2 MeV. These structures are not well described 
by any of the models used for normalization. They might be 
due to shell effects and/or a-clustering effects (since 44 Ti is 
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FIG. 5: (Color online). Comparison of experimental primary y-ray spectra (filled squares) and the ones obtained from multiplying the extracted 
p and & functions (red line) for several excitation energies. The experimental and calculated spectra are shown for excitation-energy bins of 
467 keV. 



an N = Z nucleus). An increase in the level density could 
also indicate the breaking of a nucleon Cooper pair and/or the 
crossing of a shell gap. 

We note that both the ground state and the + state at E = 
1.9 MeV are not very pronounced in Fig. [6] Also, we see that 
there is less direct decay to the ground state and the excited 
+ state compared to the 2 + state at 1083 keV and the 4 + 
state at 2454 keV, especially at higher excitation energies (see 
Fig. |3J. This could imply that there are relatively few spin- 
1 states populated in the (p,t) reaction (assuming that dipole 
radiation is dominant in this region). This is not surprising 
because the spin distribution is expected to have a maximum 
for J = 3 - 4. 

We see that the resolution on the level density is rather poor 
at low excitation energies, roughly 600 keV for the first ex- 
cited 2 + state. This may be explained by the fact that the 
level density in this region is mainly determined by decay in- 
volving high-energy y rays, which have a resolution of up to 
fa 300 keV (similar to the triton-energy resolution). As a con- 
sequence, the resolution of the level density gets better and 
better for increasing excitation energy. In addition, the rather 
large amount of other Ti isotopes in the target may give a 
smoothing effect on the extracted quantities (both the NLD 
and RSF). 



B. Radiative strength function 

Because there are no experimental data on neither the level 
spacing D nor the total, average radiative width for s-wave 
resonances (Tyo) for 44 Ti, we must find the absolute normal- 
ization parameter B of the RSF by other means. In order to 
have a rough approximation of the absolute strength, we have 
looked at experimental values of \Tyo) for other Ti isotopes 
found in Ref. Il26l . see Fig. [7] From these values, the educated 
guess of (r^o) = 1200(600) meV seems to be reasonable for 
44 Ti. The RSFs of 44 Ti for the two choices of level-density 
normalization are shown in Fig. [8] Also the upper and lower 
limits are indicated, corresponding to (Tyo) = 1800 meV for 
the CT normalization and (r^o) = 600 meV for the GHK nor- 
malization, respectively. 

From Fig. M it is seen that for E y > 4 MeV, the RSF of ^Ti 
seems to reach a relatively smooth behavior with increasing 
strength as a function of Ey. Naturally, as seen from Eq. [5] 
the slope varies depending on the normalization chosen for 
the level density. For energies below Ey ps 4 MeV, we observe 
on average a slight increase in strength for decreasing y-ray 
energy. However, we see that the data in this region display 
quite large variations, which could be due to Porter- Thomas 
fluctuations |[T9l . 
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FIG. 6: (Color online) Experimental data points (filled squares) 
normalized to the CT model (dashed line) and data points (blue, 
open circles) normalized to the calculation of GHK (dashed-dotted 
line). The black, solid line represents the known levels taken from 
Ref. I25l . 



FIG. 7: Experimental average radiative width (filled squares) for 
47_51 Ti as a function of neutron separation energy, and estimated 
width for 44 Ti (open square). 



A low-energy increase in the RSF data has been seen previ- 
ously in several light and medium-mass nuclei with the Oslo 
method E7T43T1 . with the two-step cascade method following 
neutron capture [27], and recently also in proton capture [32]. 
However, for 44 Ti, the increase at low energies is not as strong 
as in. e.g.. 56 - 57 Fe l27l. 

As of today, it is not clear whether the low-energy increase 
is due to some sort of collective decay mode(s) or if it is an 
effect of other structural effects in these nuclei. An analysis 
of simulated data using the DICEBOX code [33 1 has demon- 
strated that for light nuclei, the spin distribution of the ini- 
tial populated levels may have a considerable influence on the 
possible decay paths from these levels H141 . This seems to be 
due to a combination of several factors: (/) the low level den- 
sity in light nuclei at low excitation energy, (ii) restrictions 
on the possible populated spins of the initial levels, (Hi) the 
dominance of dipole radiation from highly excited levels, and 
(z'v) a rather large asymmetric parity distribution up to rather 
high excitation energies. As is shown in Ref. [ 14 1, these fac- 
tors may lead to a significant increase in the extracted RSF for 
low y-ray energies compared to the input RSF function used 
to generate the y-ray spectra. It is not unlikely that a similar 
effect may be present in 44 Ti as well. 



C. Models for the radiative strength function 

In general, the y decay in quasicontinuum is expected to be 
dominated by electric dipole radiations. Also, there is experi- 



mental evidence that a giant magnetic dipole resonance (also 
known as the magnetic spin-flip resonance) is present in sev- 
eral light and medium-mass nuclei (see, e.g., Ref. |34|). 

We have applied the commonly used Generalized 
Lorentzian (GLO) expression ||35l [3lH for the El strength. 
This model is given by 11361 

/GLo(£y, Tf) = - } 2 a E iT E \ x (7) 

E y r{E 7 J f ) r(g y = 0,7» " 

{E^-El.f+E^EyJff ■ El, y 

where the Lorentzian parameters Ye\, Ee\ and Gei corre- 
spond to the width, centroid energy, and peak cross section of 
the giant electric dipole resonance (GDR). We have made use 
of the parameterization of RIPL-2 ll26l to estimate the GDR 
parameters as these are unknown experimentally. Due to the 
dynamic ground-state deformation of 44 Ti (fa = 0.27 l26ll ), 
the GDR is assumed to be split in two components and thus 
two sets of GDR parameters are applied (see Tab. [I]). In addi- 
tion, we have assumed a constant temperature 7V of the final 
states in accordance with the Brink hypothesis iflTl . 

In order to account for the small increase in strength at low 
energies, we have added a Lorentzian resonance of the form 

1 ^up^y^up 

/up(£r)= 3^V (£2-£2 p )2+£ y 2 lV (8) 

with centroid energy E up = 1.9 MeV, width T up = 1.3 MeV, 
and a peak cross section (7 up that will vary according to which 
NLD model that has been used for normalization. Such a 
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FIG. 8: (Color online). Normalized radiative strength functions for 
the two normalizations of the level density. The dashed line is the 
upper limit for the CT normalization, while the dashed-dotted line is 
the lower limit for the GHK normalization. 



shape of the low-energy increase seems to be in accordance 
with data and model descriptions in Refs. |37 , 38 1 for the Mo 
nuclei, and was used also in Ref. ||39l . The constant tempera- 
ture Tf is slightly varied for the two models to give the best fit 
to the data. All parameters used are given in Tab. [j] 

The lower and upper limit of the absolute normalization 
will necessarily give slightly different parameters to get the 
best fit to the data. These parameters are also given in Tab. [I] 

Another frequently used model for E 1 strength is the stan- 
dard Lorentzian (the Brink- Axel model, see Ref. |26| and ref- 
erences therein). The expression reads 



/sLo(£y) = 



i 



(TeiEvT 2 



3% 2 h 2 c 2 (E 2 -E 



E\ ■ 



El 

F 2 r 2 ' 



(9) 



This model is independent of excitation energy, in accordance 
with our assumption behind Eq. ([IJ. However, as described in 
Ref. l26l . this model is known to generally overestimate the 
value of (Tyo) and neutron-capture cross sections. 

For the magnetic transitions, a standard Lorentzian as rec- 
ommended by the RIPL-2 library 11261 and shown in Fig.[9]is 
adopted (see Table[j]for the corresponding parameters). 

The resulting models and the extracted data for the two 
adopted NLD normalizations are shown in Fig. [9] (the mod- 
els corresponding to the best fit for the lower and upper nor- 
malization limits are not shown). We see from Fig. [9] that 
the shape of the SLO model does not fit the extracted RSF 
very well; in particular the low-energy part is very different in 
shape. The enhancement of the RSF data at low y-ray energies 
may be explained by the spin distribution of the initial levels 
as described previously. 



10" 



> 



c 
o 

1 10- 7 

c 

3 

I, 10 " 8 
C 
0) 

(0 



10- 1 



> 



c 
o 

Z 10-7 
c 

3 



- 10"' 
c 



io- ; 



10" 11 



44 Ti, CT norm 
46 Ti(y,abs) 
SLO E1 
GLO E1 
M1 
up 

GLO total 



(a) 




„Ti, GHK norm. 
46 Ti(y,abs) 
SLO E1 
GLO E1 
M1 
up 

GLO total 



(b) 




6 8 10 12 14 16 18 

y-ray energy E (MeV) 



FIG. 9: Radiative strength functions and fitted models for (a) the CT 
and (b) the GHK normalization. Note that the SLO model has not 
been fitted to the data. 



D. Temperature dependence of the strength function 

As discussed already, we rely on the Brink hypothesis when 
extracting the NLD and RSF. We would like to investigate 
if this hypothesis is reasonable, and we have therefore ex- 
tracted the RSF for two different excitation-energy regions, 



4.5 < E < 7.1 MeV and 7.3 < E < 9.9 MeV (see Fig. [TO 
for the result using the CT normalization). We observe that 
there are rather large, local fluctuations, e.g. at Ey m 5.5 MeV, 
that makes it hard to draw any firm conclusion whether the 
extracted RSF is dependent on excitation energy or not. As 
already mentioned, we expect large differences due to Porter- 
Thomas fluctuations in the decay strength of this nucleus. 
However, the gross features seem to be quite similar for the 
two excitation-energy ranges. 

To get a better understanding on the possible temperature 
dependence, we have also used the standard GLO model 
with a varying temperature corresponding to the accessible 
final excitation energies of the two ranges. The tempera- 
ture is estimated by Tf °< y / £/, and the models displayed in 
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TABLE I: Parameters used for the RSF models. 



Model 


£"£1,1 


Ofil.l 




££1,2 


Of 1.2 


T£l,2 


Tf 


Em I 


0~M1 


Tmi 




Gup 


r\ip 




(MeV) 


(mb) 


(MeV) 


(MeV) 


(mb) 


(MeV) 


(MeV) 


(MeV) 


(mb) 


(MeV) 


(MeV) 


(mb) 


(MeV) 


CT+GLO 


17.23 


43.16 


5.98 


21.58 


21.54 


9.19 


0.50 


11.6 


0.8 


4.0 


1.9 


0.060 


1.3 


CT+GLO, upper limit 


17.23 


43.16 


5.98 


21.58 


21.54 


9.19 


0.80 


11.6 


0.9 


4.0 


1.9 


0.070 


1.3 


GHK+GLO 


17.23 


43.16 


5.98 


21.58 


21.54 


9.19 


0.40 


11.6 


0.7 


4.0 


1.9 


0.015 


1.3 


GHK+GLO, lower limit 


17.23 


43.16 


5.98 


21.58 


21.54 


9.19 


0.15 


11.6 


0.3 


4.0 


1.9 


0.010 


1.3 



o 



exp., E=4.5-7.1 MeV 
exp., E=7.3-9.9 MeV 
exp., E=4.5-9.9 MeV 
std. GLO, E=4.5-7.1 MeV 
, E=7.3-9.9 MeV 
, E=4.5-9.9 MeV 
T=0 

E=9.9 MeV 



std. GLO 
std. GLO 
std. GLO; 

std. GLO 




y-ray energy E^ (MeV) 

FIG. 10: (Color online) Radiative strength functions extracted for 
two different excitation-energy regions as compared to the one from 
the total excitation-energy region under consideration (all cases are 
normalized to the CT level density). The standard GLO model is also 
shown for these ranges of excitation energy and the extreme cases 
T f = and E = 9.9 MeV. 



Fig. 10 represent the average RSF within the three excitation- 
energy regions 4.5 < E < 7.1 MeV, 7.3 < E < 9.9 MeV, and 
4.5 <E< 9.9 MeV. The maximum final temperature reached 
is in the range 0.50-1.14 MeV for initial excitation energies 
between 4.5 - 9.9 MeV and with £™ n = 2.2 MeV. 

We observe that for y energies between 2.2 — 5.2 MeV, 
the upper range can give more than a factor of two larger y 
strength in the GLO model than the lower range. The largest 
experimental fluctuations are also of the order of a factor two, 
but on average, the experimental data seem not to have such a 
strong dependence on the final excitation energy as the stan- 
dard GLO model would imply. 

Finally, we have calculated the GLO model for the direct 
decay to the ground state {Tf = MeV) and for the high- 
est temperatures reached, namely for initial excitation energy 
E =9.9 MeV. These cases are also displayed in Fig. 10 It 



is clear that the zero-temperature calculation is not similar to 
the experimental data below £|Rj6 MeV. In fact, the experi- 
mental data lie in between the two extreme cases. As argued 



above, the overall shape of the RSF seem to be rather similar 
for the two excitation-energy ranges. Thus, it is probably quite 
reasonable to apply a constant temperature in the GLO model 
in accordance with the Brink hypothesis. We note however 
that the GLO model averaged over the total excitation-energy 



range covered (blue solid line in Fig. 10 1 agrees rather well 



with the average shape of the experimental data. 



IV. REPRODUCTION OF EXPERIMENTAL PRIMARY 
y-RAY SPECTRA 

As described in the previous section, both the NLD and 
RSF deduced from the present experiment are affected by 
severe uncertainties associated with the normalization proce- 
dure that in the case of 44 Ti cannot be reliably constrained 
by additional experimental data. For this reason, the different 
NLD and RSF models are now directly tested on the primary 
y-ray spectra, which in turn correspond to the fundamental 
quantity entering the description of the radiative decay in re- 
action models. 

There are 24 experimental spectra for excitation energies 
4.5 < E < 9.9 MeV with bin size 233.4 keV. Due to the poor 
statistics, we have compared the average of two bins as in 
Fig. [5] For the level density, we have applied the known levels 
for E < 3.7 MeV, and either the CT model or the GHK calcu- 
lation above this energy. The calculated spectra are scaled to 
get the best possible agreement with the experimental spectra, 
which are normalized such that _ E mm P(E,Ey) = 1. 

We have also calculated the primary y spectra using the 
standard GLO model with a variable temperature in combi- 
nation with the two level density models. For each initial 
excitation energy (4.5 < E < 9.9 MeV), we have used the 
GLO model with Tf « \fEj and made an average of all these 
RSFs for the whole excitation-energy region that is used for 
the analysis. 

Finally, we have applied the SLO model in combination 
with the two NLD models. 



The six model combinations shown in Figs. 11-12 thus cor- 
respond to: 

• input 1 : the CT level density and the fitted GLO model 
with constant temperature. 

• input 2: the GHK level density and the fitted GLO 
model with constant temperature. 

• input 3: the CT level density and standard GLO model 
averaged over the experimental excitation-energy range 



(blue, solid line in Fig. 10 1 
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input 4: GHK level density and standard GLO model 
averaged over the experimental excitation-energy range 
(blue, solid line in Fig.fTO). 



• input 5: the CT level density and the SLO model. 

• input 6: the GHK level density and the SLO model. 

We find that for the lower excitation energies (E < 7.0 
MeV), all the input models give a rather good reproduction of 
the spectra, and also they give relatively similar results. How- 
ever, for the higher excitation energies (7.4 < E < 9.8 MeV), 
there are clear deviations for the different inputs, and input 
4 gives a significantly worse fit than the others. This is not 
so surprising, considering that the slope of the standard GLO 
model is very different from the RSF model fitted to our data 



using the GHK level density. It is seen in Fig. 12 that this mis- 
match in slope leads to a wrong overall shape of the primary 
spectra (overestimating the low-energy part and underestimat- 
ing the high-energy part). The models that are consistent in 
slope behave much better. 

In general, for the high-energy region, input 2 and input 3 
give the best reproduction of the experimental spectra out of 
the six model combinations. The SLO model gives also rather 
good results, in particular in combination with the GHK level 
density. This could be an indication that the overall shape of 
this model might be correct, although the absolute value is 
probably too large (as seen in Fig. 15). However, in combina- 
tion with the CT level density the spectra calculated with this 
model underestimate the intensity of the low-energy y rays 
compared to the measured ones. 

Naturally, using the extracted p and £? data points give a 
much better fit than any model, as structures especially in the 
level density are lost when using the smooth models. 



V. CAPTURE CROSS SECTION AND 
MAXWELLIAN-AVERAGED RATE 

The model combinations described in Sec. |IV] can fur- 
ther be applied for estimating the 40 Ca(a, y) 44 Ti cross sec- 
tion and the corresponding reaction rate. We have used the 
code TALYS [24 1 for the cross-section and reaction-rate cal- 
culations. In all calculations we have used the a optical- 
model potential (OMP) of McFadden and Satchler g(j]. We 
have also tested other a-OMPs such as the one developed by 
Demetriou, Grama and Goriely [41 1; however, it turns out that 
in this case, the results are rather insensitive to the choice of 
a-OMP. 

The model sets applied for the cross-section and reaction- 
rate calculations are: 

• input 1 : the CT level density and the corresponding fit- 
ted GLO model. 

• input 2: the GHK level density and the corresponding 
fitted GLO model. 

• input 5: the CT level density and the SLO model. 

• input 6: the GHK level density and the SLO model. 



• input 7: the CT level density and the standard GLO 
model. 

• input 8: the GHK level density and the standard GLO 
model. 

• input 9: the CT level density and the corresponding 
fitted GLO model for the upper normalization limit 
((Lyo) = 1800 meV). 

• input 10: the GHK level density and the correspond- 
ing fitted GLO model for the lower normalization limit 
((r^o) = 600 meV). 

For the standard GLO model, we have followed the pre- 
scription used in TALYS, where the temperature of the final 
states is calculated from the relation [26| 



T f = J(E,-A-E r )/c 



(10) 



Here, £; is the initial excitation energy in the compound nu- 
cleus, A is a pairing correction, and a is the level density pa- 
rameter. There is thus a specific GLO RSF for each initial 
excitation energy (no averaging over a large excitation-energy 



window as was done with inputs 3 and 4 in Sec. IV i. 

The Maxwellian-averaged astrophysical reaction rate 
Na (ctv) (T) for a temperature T is given by 



N A (av) (T) = 



1/2 n a 

Jim I (kT) 3 / 2 G(T) X 

(E + E?) 



(ID 



o y (27° 11 1 



1) 



kT 



dE, 



where m is the reduced mass of the initial system of projectile 
(here the a particle) and target nucleus (here 40 Ca), k is the 
Boltzmann constant, Na is Avogadro's number, E is the rel- 
ative energy of the target and projectile, I^,Ex are the spin 
and excitation energy for the excited states labeled ju, and 
is the reaction cross section. Further, the temperature- 
dependent, normalized partition function G(T) reads 



G(T 



l)/(2/°+l)exp(-/47fcr). (12) 



See the TALYS documentation for more details [24). 

For the above expressions, local thermodynamic equilib- 
rium is assumed in the astrophysical environment, so that the 
energies of both the targets and projectiles as well as their rel- 
ative energies E obey Maxwell-Boltzmann distributions cor- 
responding to the temperature T at that location. Also the rel- 
ative populations of the various nuclear levels with spin and 
excitation energy 1^ ,E% obey a Maxwell-Boltzmann distribu- 
tion. 



A. Uncertainties in the calculations 

The key ingredients in the rate calculations (NLD, RSF, 
and a-OMP) all contribute to the uncertainties but in differ- 
ent temperature regions. By varying the a-OMP, the rate will 
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FIG. 1 1 : (Color online) Experimental primary y spectra (filled squares) and those obtained from multiplying the extracted p and ,9~ functions 
(red line) for 4.6 < E < 7.0 MeV. These are compared with calculated spectra using the six inputs as described in the text. The experimental 
and calculated spectra are given for excitation-energy bins of 466.8 keV. 



change most at low temperatures (below w 1.5 • 10 K), while 
the NLD and RSF both have a significant impact on the rate 
at higher temperatures (above « 1.5 • 10 9 K). 

The 40 Ca(a,y) reaction only populates states with total 
isospin zero in 44 Ti, and dipole transitions with no change 
in total isospin are suppressed in self-conjugate N ~ Z nu- 
clei ll42l . Since complete isospin mixing is assumed in the 
determination of the RSF, a corrective factor must be included 
before applying the extracted RSFs in the cross section and re- 
action rate calculations. For a-captures, a standard prescrip- 
tion is to divide the RSF by a constant factor /; s0 of typically 
5 to 8 Ella, giving a reduced RSF /' off = f/f iso . 

However, it is not obvious how the value of this correc- 
tion factor should be determined, as it is related to the de- 
gree of isospin mixing (for a large degree of mixing the cor- 
rection factor should be small and vice versa). Also, high- 
energy 7 rays decay to low-lying states, where the isospin 
mixing might be small, while low-energy 7 rays decay in the 
quasi-continuum, where one expects a large degree of mix- 
ing. Therefore, one could in principle expect that /; so would 
vary as a function of 7 energy. To estimate such a function 
is however a very complicated task and beyond the scope of 
the present work. We will therefore assume a constant /; so , 
although this is probably quite crude. 



Another source of uncertainty is the value of (Tyo). Since 
this value is not known experimentally, the uncertainty in this 
quantity is correlated to the uncertainty in /; so . This is because 
one can obtain basically the same cross section and reaction 
rate for a range of (r^o) values by adjusting /; so correspond- 
ingly (a small value of (Tyo) in combination with a small /j so 
and vice versa). 

We have calculated the integral cross section of the 
(a, 7) 44 Ti reaction for incoming a energies between 2.1— 4.2 
MeV. This corresponds to excitation energies in 44 Ti in the 
range E = 7.2 — 9.3 MeV, which is the most relevant region 
for astrophysics. In the calculations, we have adopted a con- 
stant /; so = 5, and we have also tested a larger correction fac- 
tor of /i so = 8. In addition, we have considered the assumed 
uncertainty of 50% in the estimated value of (Tyo) combined 
with the uncertainty in slope (either the CT or the GHK level 
density, which represent the extremes in slope of the extracted 
RSF). The results for all the considered inputs are given in 
Tab. mi 

We see from Tab. that changing the isospin correction 
factor leads to a change in the calculated cross section of typi- 
cally 2-3 fib, except for the SLO model where different values 
of /i so give up to w 5 jJ.b change in (7 ave - The uncertainty in 
(r^o) and slope will also give a change in a ave of up to « 5 



11 



0.16 
0.14 
c 0.12 

§ °- 1 

| 0.08 
12 0.06 
tt 0.04 
0.02 



0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 



E = 7.4 MeV 




■ data 
— fit 
a input 1 
o input 2 
input 3 
* input 4 
□ input 5 
! nput 6 



E = 7.9 MeV 




E = 8.4 MeV 




E = 8.8 MeV 



E = 9.3 MeV 



E = 9.8 MeV 




i .... i .... i . 





2 3 4 5 6 7 8 
E y (MeV) 



2 3 4 5 6 7 
E Y (MeV) 



2 3 4 5 6 7 8 
E Y (MeV) 



FIG. 12: (Color online) Same as Fig.[TT]for 7.4 < E < 9.8 MeV. 



TABLE II: Integral cross sections for the various model combina- 
tions. 



Model combination a ave (M^) 

/iso 5 _/[ so 8 


input 1 


7.1 


4.9 


input 2 


6.9 


4.9 


input 5 


17.8 


13.2 


input 6 


19.2 


14.5 


input 7 


5.6 


3.8 


input 8 


6.9 


4.9 


input 9 


9.7 


6.9 


input 10 


5.0 


3.4 



Mb. 

We therefore conclude that the absolute value of the reac- 
tion cross section of the a capture on 40 Ca is highly uncer- 
tain. Also, intrinsically very different models of the level den- 
sity and the RSF may yield practically the same cross section 
by adjusting /; so and/or (r^o). Thus, we find that although 
our data may put a constraint on the functional form of the 
NLD and RSF, and in particular the correlated slope of the two 
quantities, other data are needed in order to further constrain 
the cross section and reaction rate. This will be addressed in 
the following. 



B. Comparison with other data 

We have compared our calculations with two recent data 
sets, one from Nassar et al. |6l and one from Vockenhuber 
et al. H). In the following discussion we have used /; so = 5, 
unless stated otherwise. 

Nassar et al. |6| performed an integral measurement on the 
40 Ca(a, 7) 44 Ti cross section corresponding to an energy win- 
dow of E a —2.1 —4.2 MeV for the incoming a particles. 
The originally estimated energy-averaged cross section was 
Cavi? = 8.0(11) fib; however, it was pointed out by Vocken- 
huber et al. il l that due to an overestimate of the 40 Ca stop- 
ping power, the cross section should be reduced by about 10%: 
cC P ~ 7.0(10) Mb. 

Using the model combinations input 1 and input 2, we get 
Cave =7.1 and 6.9 Mb for the CT and GHK normalization, 
respectively (see Tab. In}, This is in excellent agreement with 
the Nassar results. Note that if we exclude the low-energy en- 
hancement in input 1 and 2, (J a ve yields 6.6 and 6.7 Mb for 
the CT and GHK normalization, respectively. Thus, the ef- 
fect of the small low-energy enhancement is not significant 
compared to neither the current experimental uncertainty on 
the integral cross section, nor the uncertainties in the calcula- 
tions due to the isospin suppression factor and the unknown 
absolute normalization of the RSF. 

Further, if we apply the standard GLO model in the calcula- 
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FIG. 13: (Color online) (a) Data on the Maxwellian-averaged reac- 
tion rate for the 40 Ca(a,y) reaction from Ref. (T| are compared to 
calculations with inputs 1, 2, 7 and 8. (b) Ratio of the calculated 
reaction rates and experimental data. The upper dotted line indicates 
a 40% deviation from the experimental data, while the lower dotted 
line indicates a 20% deviation. 



tions, we obtain rj aV e =5.6 and 6.9 /zb using the CT and GHK 
level density (input 7 and input 8), respectively. The former 
is slightly smaller than the lower experimental limit on ofve ■ 
The latter agrees very well with the Nassar data; however, we 



noted in Sec. IV that the GHK level density in combination 
with the averaged, standard GLO did not reproduce our ex- 
perimental spectra above £w7 MeV. Therefore, input 8 is 
probably not correct. 

The calculations with the SLO model are many standard 
deviations too large compared to the Nassar data, even for 
/i so = 8. Thus, we will not use inputs 5 and 6 further. We 
also see that the upper normalization limit, input 9, is only 
acceptable for /i so = 8. On the other hand, the lower normal- 
ization limit (input 10) seems to be a bit too low for /; so = 5 
and even more so for /; so = 8. 

With the inputs 1,2,7 and 8, we have estimated the reaction 
rate according to Eq. (10), and compared to the DRAGON 
data of Vockenhuber et al. (TJ. The results are shown in the 
upper panel of Fig. [13] As seen here, all these four inputs give 
a reasonable reproduction of the measured rates, and they can 
hardly be distinguished from each other over the entire tem- 



perature range. We also note that the cross section range spans 
over ten orders of magnitude, and it is quite impressive how 
well the calculated rates follow the overall functional form of 
the data. 

However, we observe that some models overestimate the 
rate somewhat, especially at higher temperatures. Now, the 
experimental rate could be too low for temperatures above 
Tg 4 MeV because of missing resonances (only six reso- 
nances with measured resonance strengths are known for ex- 
citation energies at and above 9.3 MeV HI). We note however 
that the rate at these high temperatures is not important for the 
final 44 Ti mass yield Q). 

The combination of CT level density and standard GLO (in- 
put 7) gives a better reproduction of the DRAGON data at high 
temperatures than the others. However, this combination gave 
a too small average cross section as compared to the Nassar 
result. We therefore confirm the slight inconsistency between 
the two measurements, as discussed in Ref. [1|. As already 
mentioned, we observed that the combination of the GHK 
level density and the averaged, standard GLO model (input 
4) was not able to give a reasonable description of our pri- 
mary 7 spectra for initial excitation energies above s» 7 MeV. 
Therefore, input 8 cannot be recommended on the basis of the 
present data. 

In order to get a clearer picture of the deviation of the calcu- 
lated rates vs. the DRAGON data, the ratio of the calculated 
and experimental rate is shown in the lower panel of Fig. 13 
Here, it is seen that input 7 gives the best fit to the DRAGON 
data for all temperatures above Tg m 1.7 (Tg = 10 9 K). The 
large overestimate of the rate for Tg < 1.5 is common for all 
model predictions and is due to problems with the a OMP at 
very low energies. 



It is seen from the lower panel in Fig. 13 that all the calcu- 
lated rates lie within 40% of the DRAGON data for 1 .8 < Tg < 
4.5, and for the 20% upper limit in the range 2.2 < Tg < 3.7. 
However, only input 7 (CT level density and standard GLO) 
is within the experimental error bars for 1.8 < Tg < 4.5. 

It should be noted that by using the lower normalization 
limit for the RSF (blue, dashed-dotted line in Fig. [8]), an excel- 
lent agreement with the DRAGON data is obtained. However, 
as stated previously, a good reproduction of the DRAGON 
data implies a too low integral cross section (in this case 
Cave =5.0 /J.b) as compared to the Nassar data. On the other 
hand, using the upper normalization limit for the RSF (dashed 
line in Fig. [8}, it is necessary to use the larger value of /j so = 8 
to obtain a reasonable agreement with the data (see Tab.[n|. 

To summarize, we have found that the model combina- 
tions input 1, input 2, and input 8 are in excellent agreement 
with the Nassar cross-section measurement. However, input 8 
(GHK level density and standard GLO) is not in accordance 
with our primary y spectra for the relevant excitation-energy 
region (E = 7.2 — 9.3 MeV) and should therefore not be used. 
The DRAGON data are best described with input 7 (CT level 
density and standard GLO) and the lower normalization limit 
of our RSF data (input 10). We also see that the small en- 
hancement in the RSF at low Ey is not very important for the 
integrated cross section or the rate, as it gives a contribution 
of maximum 0.5 /J.b. 
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VI. SUMMARY 

Particle-y coincidence data of the 46 Ti(p,f y) 44 Ti reaction 
have been measured at OCL. By use of the Oslo method, pri- 
mary y-ray spectra have been extracted for initial excitation 
energies in the range E = 4.5 — 9.9 MeV. From these spectra, 
the functional form of the level density and the RSF have been 
determined. 

We have shown that a consistent normalization of the NLD 
and RSF is necessary in order to obtain a reasonable repro- 
duction of the primary spectra. Also, the RSF seems to be 
independent of excitation energy and thus of the tempera- 
ture, in accordance with the Brink hypothesis. However, the 
GLO model with a variable temperature, and averaged over 
the excitation-energy region, gives a rather good description 
of the overall shape of our RSF data when normalized to the 
CT level density. 

Using input models consistent with our data give an excel- 
lent reproduction of the Nassar integral measurement of the 



Ca(a,y) reaction cross section. Also the DRAGON data of 
the reaction rate are rather well reproduced. However, we note 
that there is a discrepancy between the two datasets, so that it 
is not possible to use one model combination to obtain optimal 
agreement for both measurements simultaneously. Neverthe- 
less, on the basis of our present data, it is clear that certain 
model combinations are not acceptable, although they might 
give reasonable results compared to the cross-section data. 
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